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Ren  Sicong  (Northwestern  Polytechnical  University)  and 
Zhao  Yongkang  (Baocheng  Instruments  Factory) 


Abstract 


This  paper  proposes  a  new  inertial  navigation  system  com¬ 
bining  the  platform  type  and  strap-down  type  inertial 
navigation  systems.  In  this  system,  the  azimuth  rate  platform 
without  azimuth  stabilized  loop,  azimuth  coordinate  resolver 
and  synchronizer  are  used.  The  azimuth  angles  of  the  platform 
and  vehicle  are  obtained  by  an  integrator  from  azimuth  rate 
signals  measured  by  an  azimuth  gyroscope  supported  on  a  hori¬ 
zontal  gimbal.  This  type  of  inertial  navigation  system  is 
suitable  for  vehicles  without  large  angle  pitching  maneuvers 
such  as  transports,  aerodynamic  and  ballistic  missiles  etc. 

This  paper  discusses  the  operational  principles  of  the  azimuth 
rate  platform,  mechanization  equations,  distinguishing  features 
of  initial  alignment,  and  the  calibration  and  compensation  for 
drifts  of  gyroscopes;  at  the  same  time,  error  propagation 
characteristics  caused  by  various  major  error  sources  for  navi¬ 
gation  positioning,  velocity  and  attitude  are  simulated  on  a 
computer.  In  the  conclusion,  we  point  out  that  the  simplicity 
of  the  platform  structure,  smail  volume  and  weight,  high  relia¬ 
bility,  and  the  possibility  of  calibration  and  compensation  for 
drifts  of  the  azimuth  gyroscope  are  the  outstanding  advantages 
of  this  inertial  navigation  system.  Furthermore,  if  we  use  a 
special  optical  system  coordinated  with  known  azimuth  angles 
and  latitudes  of  land  marks,  it  is  not  only  possible  to  realize 
fast  alignment  but  also  the  calibration  and  compensation  of 
horizontal  gyroscopes. 


I. 


Putting  Forward  the  Problem  and  Brief  Introduction  of  the 
Platform  Structure 


1 


Even  though  there  are  presently  many  types  of  inertial 
navigation  and  guidance  systems  widely  used  in  aeronautics, 
marine  navigation  as  well  as  aerodynamic  guided  missiles,  yet 
they  can  basically  be  summed  up  in  the  two  large  categories  of 
the  platform  type  and  strap-down  type. 

The  basic  distinguishing  features  of  the  platform  type 
system  are  that  it  uses  a  stabilized  loop  to  separate  the  move¬ 
ments  of  the  inertial  device  passing  the  gimbal  and  the  vehicle 
thus  causing  them  to  be  in  an  excellent  operating  environment 
and  their  dynamic  ranges  to  be  relatively  small.  It  especially 
causes  the  influence  of  the  gravitational  acceleration  to  de¬ 
crease  to  minimum  when  the  semi-analytic  type  system  is  operating 
in  a  horizontal  coordinate  system.  At  the  same  time,  the  com¬ 
puted  relationship  is  also  relatively  simple  which  proves  the 
high  accuracy  of  the  system.  However,  the  major  drawbacks  of  the 
system  are  that  the  structure  of  the  system  is  complex,  manu¬ 
facturing  costs  are  very  high,  the  volume  and  weight  are  rela¬ 
tively  large  and  at  the  same  time,  the  number  of  gimbals  and 
slip  rings  are  quite  numerous  so  that  the  reliability  is  somewhat 
affected  etc. 

The  inertial  device  of  the  strap-down  system  is  directly 
strapped  on  the  vehicle  and  thus  its  major  advantages  are  its 
structure  is  simple,  reliability  good  and  the  attitude  information 
which  does  not  require  an  electromechanical  component  for  switch¬ 
ing  over  has  relatively  high  precision;  however,  the  operating 
environment  of  the  inertial  device  is  adverse  especially  when  the 
vehicle  carries  out  sharp  maneuvers  as  well  as  for  their  dynamic 
range  which  must  be  very  wide  and  the  extremely  harsh  require¬ 
ments  for  the  gyroscope.  Furthermore,  because  the  direction  of 
the  inertial  device  facing  the  gravitational  field  is  constantly 
changing,  the  compensation  technique  for  its  errors  and  gravita¬ 
tional  force  is  quite  complex;  requirements  for  the  word  length, 
speed  and  storage  capacity  of  the  computer  are  relatively  high. 
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Therefore,  we  encounter  very  large  difficulties  in  guaranteeing 
the  operating  precision  of  the  strap-down  system. 

For  a  large  number  of  vehicles,  because  of  structural, 
mechanical  and  other  specific  reasons,  when  maneuvering,  the 
angular  rate  difference  of  the  roll  and  yaw  movements  is  very 
far.  For  example,  for  common  aircraft,  when  the  roll  angular 
rate  can  reach  or  exceed  360  degrees  per  second,  the  yaw  angular 
rate  is  commonly  only  2-3  degrees  per  second.  The  maximum  cannot 
exceed  10  degrees  per  second.  In  this  way,  if  we  say  we  want  to 
use  a  conventional  inertial  gyroscope  to  accurately  measure  the 
maximum  angular  rate  of  10  degrees  per  second  and  not  be  able  to 
have  very  great  difficulties,  it  will  be  very  difficult  to  use 
it  to  measure  angular  rates  above  360  degrees  per  second.  For 
just  this  reason,  in  order  to  fully  bring  into  play  the  advan¬ 
tageous  points  of  the  platform  and  strap-down  inertial  naviga¬ 
tion  system,  we  propose  a  new  azimuth  rate  inertial  navigation 
system  which  can  simplify  the  platform,  lower  the  manufacturing 
costs  and  raise  reliability.  The  platform  which  it  uses  elimin¬ 
ated  the  azimuth  stabilized  loop,  azimuth  coordinate  resolver 
and  azimuth  synchronizer  of  the  common  three  axis  platform.  The 
azimuth  angles  of  the  platform  and  vehicle  are  derived  by  an 
integrator  from  azimuth  rate  signals  measured  by  an  azimuth  rate 
gyroscope  supported  on  a  horizontal  gimbal. 

When  this  type  of  system  was  first  brought  forth  in  1978, 
it  was' called  the  "strap-down  azimuth  inertial  navigation  system" 
or  the  "semi-strap-down  type  inertial  navigation  system." 
Generally  speaking,  in  order  to  simplify  the  platform,  we  can 
also  have  a  "strap-down  pitching  system"  and  even  an  "azimuth 
and  pitching  strap-down  system"  which  only  has  a  roll  stable 
axis.  Howevdr,  when  we  consider  applying  this  system,  it  is 
most  advantageous  in  a  horizontal  coordinate  semi-analytical 
inertial  navigation  system  with  a  not  very  large  azimuth  angle 
rate;  it  is  basically  the  same  as  the  "rotating  azimuth*, 


"free  azimuth",  "drifting  azimuth”  and  other  systems.  Further, 
when  the  platform  has  attitude  errors,  the  azimuth  gyroscope  is 
also  not  completely  "strapped"  with  the  vehicle.  In  order  to 
lose  the  meaning  of  " strapped -down"  we  therefore  further  changed 
the  name  to  "azimuth  rate"  or  "analytical  azimuth"  inertial 
navigation  system. 


The  configuration  of  the  azimuth  rate  platform  is  shown  in 
Fig.  1.  It  can  use  two  dual  degrees  of  freedom  gyroscopes  or 
one  single  degree  of  freedom  gyroscope.  It  can  also  be  composed 
of  one  single  degree  of  freedom  angular  rate  and  two  single 
degree  of  freedom  integrating  gyroscopes. 


Fig.  1  Configuration  of  the  azimuth  rate  platform. 


\ 


Fig.  2  Scheme  of  the  azimuth  rate  platform. 

Fig.  2  shows  the  principle  of  the  azimuth  rate  platform 
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composed  from  two  dual  degrees  of  freedom  gyroscopes.  We  can 
see  from  the  figure  that  the  two  dual  degrees  of  freedom  gyro¬ 
scopes  2  and  6  and  the  three  accelerometers  10,  11  and  15 
installed  on  platform  9  form  the  inertial  measuring  component 
supported  by  horizontal  gimbal  17  and  that  the  horizontal  gimbal 
is  supported  by  the  stand  strapped  together  with  vehicle  1.  The 
main  axis  of  gyroscope  2  is  in  the  plumb  state  during  normal 
operations  and  it  uses  corresponding  amplifiers  3  and  14  as  well 
as  stabilized  electric  machines  16  and  13  to  realize  stabiliza¬ 
tion  of  the  platform's  pitch  axis  and  roll  axis;  the  main  axis 
of  gyroscope  6  crosses  with  gyroscope  2  and  uses  amplifiers  5 
and  8  to  realize  feedback  self-locking  and  operates  in  a  rate 
condition.  Because  the  feedback  current  of  the  passage  surround¬ 
ing  the  azimuth  axis  is  the  measurement  of  the  azimuth  angle 
changing  rate,  integrator  4  can  compute  the  azimuth  angle. 

Because  the  supporting  axis  line  of  the  horizontal  gimbal  and 
the  longitudinal  axis  of  the  vehicle  are  parallel  to  eachother, 
angular  transucer  7  can  put  out  roll  attitude  signals;  cor¬ 
respondingly,  we  obtain  the  pitch  attitude  signals  from  angular 
transducer  12.  It  is  very  clear  that  this  type  of  structure 
simplifies  the  platform  into  a  simple  device  which  aside  from  the 
azimuth  gyroscope  only  has  a  horizontal  gimbal  similar  to  that  of 
the  vertical  gyroscope.  However,  when  in  pitch  angles  larger 
than  60  degrees,  the  maneuvering  movements  of  the  vehicle  pos¬ 
sibly  influence  the  normal  operation  of  the  platform.  In  order 
to  adapt  to  this  type  of  situation,  it  is  necessary  to  increase 
the  external  rolling  ring  and  corresponding  stabilizing  loop. 

If  the  azimuth  rate  inertial  navigation  system  is  not  used 
for  vehicles  which  make  complete  degree  of  freedom  rolling  and 
pitching  movements,  the  entire  platform  need  not  use  the  conduct¬ 
ing  slip  ring  but  completely  use  the  soft  lead  wire.  This  not 
only  simplifies  the  structure  of  the  platform  and  makes  main¬ 
tenance  convenient  but  even  more  important  it  can  also  raise 
reliability. 
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Because  the  operating  state  of  the  azimuth  rate  gyroscope 
is  similar  to  that  of  the  strap-down  system,  when  compared  with 
the  azimuth  gyroscope  of  the  common  platform,  its  calibration 
system  errors  and  nonlinearity  are  possibly  enlarged  because 
the  dynamic  range  is  relatively  large;  however,  because  the 
vibration  isolation  of  the  horizontal  gimbal  and  the  platform  are 
basically  in  horizontal  states,  the  gravitational  effects  are 
relatively  small  and  it  is  easy  to  guarantee  the  operating  pre¬ 
cision  of  the  azimuth  gyroscope.  When  we  use  the  pulse  rebal¬ 
ance  loop,  this  can  also  avoid  the  additional  switching  errors. 

II.  Mechanization  Equations  and  Structural  31ock  Diagram 

We  select  the  north,  west  and  vector  geographical  coordinate 
system  NWV.  As  shown  in  Fig.  3,  ideal  platform  system  X^,Y_Z_ 
follows  the  vehcile  rotating  around  vertical  axis  V  and  the 
included  angle  of  the  X^  relative  meridian  N  is  platform  azimuth 
angle  (V . 


Fig.  3  Coordinate  system. 

Based  on  the  specific  power  equation,  we  determine  the 
velocity  of  the  platform  relative  to  the  earth  from  the  following 
differential  equations: 


V  *  -  /x  -  (2Q»  -f  <By)  Vz  4-  (2 Qz*  -r  CDZ»)  Fy 
Vy™fr~  (2Qz»-rOi)z*)f/’ x  +  (2Q*  -f-<ax)^ z 
V>“/*-(2Qx  J-ffi1)^r+(2QT  +  aiY)^x-  <7 


In  the  equations,  is  the  specific  power  measured  by  c 
corresponding  axis  accelerometer;  S2 .  is  the  earth's  c  oon- 
ent;  60  is  the  angular  velocity  component  of  the  platf  n 
relative  to  the  earth's  rotation;  g  is  the  gravitation* 
acceleration;  and  symbol  b  respresents  the  angular  velocity  of 
the  platform  driven  by  the  vehicle. 


In  order  to  carry  out  navigation  computations,  we  can  use 
the  relationship  of  the  directional  cosine  matrix  (C)  formed 
from  longitude  A  /  latitude  L  and  azimuth  angle  cu  between  the 
platform  and  earth  system  and  its  change  rate,  and  the  (x>. 

Cc3-CCK«0  (2) 

In  the  formula 


'  sin  X  sii 
—  T 


—sin  I  cos  X  cos  a 
—cos  X  sin  a 

—sin  I  sin  X  cos  a 
cos  L  cos  a 


sin  X  cos  a 

sin  I  cos  X  sin  a 

—  cos  X  cos  a 

sin  L  sin  X  sin  a 

—  cos  L  sin  a 


cos  L  cos  X 
cos  £  sin  X 


It  is  very  clear  that  bD  and  bb>  are  derived  based  on  V 

Ax  Y 

and  V..  and  the  earth's  curvature  is  computed  by  computer  and 


bringing  pressure  on  the  gyroscope  causes  the  platform  to  main¬ 
tain  a  control  quantity  of  a  horizontal  state.  Thus,  £D7, 

Li 

requires  that  ^e^-ete<^  from  azimuth  rate  gyroscope 

feedback  control  angular  rate  t^iat 

Further,  when  we  consider  the  independent  mixed  attitude 
system  formed  from  the  atmospheric  pressure  altitude  information 
and  vertical  accelerometer  information,  we  can  show  the  mechan¬ 
ization  relationship  of  the  azimuth  rate  inertial  navigation 
system  in  the  block  diagram  of  Fig-.  4. 


A. 


Fig.  4  Block  diagram  of  the  azimuth  rate  inertial  navigation 
system. 

and  Gj_  are  the  accelerometer  and  gyroscope  of  the 
corresponding  axis;  ha  is  the  etmospheric  pressure 
altitude  information:  and  are  the  pitch  and 

roll  state  information;  /l0  and  L0  are  the  initial 
longitude  and  latitude;  Q  is  the  rotation  angle  rate 
of  the  earth. 


III.  Special  Features  of  the  Initial  Alignment  of  the  Azimuth 
Rate  Platform 

If  we  do  not  consider  the  structural  errors,  the  azimuth  of 
the  ideal  platform  system  is  also  the  parking  position.  For 
just  this  reason,  when  we  refer  to  the  runway's  middle  line  or 


the  special  aircraft  parking  markers  placed  on  the  aircraft 
parking  area  to  determine  the  parking  position  of  the  aircraft, 
we  not  only  possibly  realize  fast  speed  alignment  but  can  also 
carry  out  calibration  and  compensation  of  the  drifts  of  the 
gyroscopes.  In  order  to  conveniently  explain  the  problem,  we 
must  first  establish  the  attitude  error  equation.  It  is  consid¬ 
ered  that  when  the  attitude  error  angle  is  small  in  number  we 
can  obtain 

^  =  a,-C£3  r  (6) 

In  the  equation,  CP  is  the  attitude  error  angle  vector; 

43  T=  £?=[£i?x  £j£?Y  is  the  angular  rate  vector  of  the  ideal 

platform  system; 

'  1  0  -CPr  <| 

Cf-  0  1  <Px  | 

*  <Pt  -«P«  1  \ 

is  the  attitude  error  matrix;  (*>  =  CO  +  43 „„+E  is  the  angular 

P  c  ZPP 

rate  vector  of  the  platform  system's  rotation  and 
—  T 

£Oc=[  43Xc  -^Yc  is  the  contro1  added  to  the  platform; 

<SDzb[o  o  S?zbp]T=[o  o  <PY  <h?x-  ^x  ^Y+  ^  is  the  an<?ular 

rate  of  the  vehicle  with  platform  rotation;  E=[E  E  0]T  is  the 
shift  angular  rate  of  the  platform. 

By  substituting  the  corresponding  relationship  into  formula 
(6)  we  obtain 

QcosLcosa  -f-'FyQsin  L  +£"* 

(pY=sQjy,  — cpxQgin  £  +  Qcxsisin  a +£y  (7) 

0  =*Qz*r  —  spyQcosLcosa  -i-  ®,Qcos  L  sin  a  —  Qsin  £ 

This  formula  shows  that  when  the  azimuth  rate  platform  is  in 
initial  alignment,  the  azimuth  attitude  errors  do  not  produce 
any  changes.  It  also  clearly  reflects  the  basic  special  features 
of  this  system's  initial  alignment.  We  can  display  formula  (7) 
as  the  block  diagram  shown  in  Fig.  5.  In  it,  coefficients 


and  directly  determine  the  dynamic  characteristics  of 
the  leveling  process.  We  shouxa  rationally  make  a  selection 
based  on  the  transient  time  and  the  basic  requirements  of  anti¬ 
interference  . 

Under  stable  conditions,  from  formula  (7)  we  obtain 
cdx,(oo)  -  Ex  —  Q cos  L  cos  a 

(8) 

=*£T+  Qcosisina 


From  this,  we  can  find  the  azimuth  angle 


—  CDt.(co 


Q cos /.sir  a  +  Er 


*  Qcos  Z.  cos  a  — ...  _ 


It  is  very  clear  that  the  platform  drift  has  direct  influence  on 
the  calculation  precision  of  the  azimuth  angle.  Assuming 
a  =a+Aa,  and  letting  <£>  (oo)  =  £?cosLcosa  and  =  QcosLsina  , 

C  XC  iC  c 

from  formula  (8)  we  can  obtain 
QcosLsinaAa-sE, 


Q  cos  L  cosa  Aa  =  Et 


If  we  separately  multiply  this  formula  by  sina  and  cosa  and  add 
them  to  eachother,  we  then  find  the  calculation  error  of  the 
azimuth  angle  is 

Exsitt  a  +ErCOga  _  £»  /■,,) 


Qcos  L 


Qcos  L 


In  the  formula,  E.  =Evsina+Evcosa  is  the  equivalent  western  drift 

W  a  Y 

of  the  platform.  This  conclusion  is  completely  identical  to 
other  systems. 


In  order  to  calculate  the  compensation  quantity  of  the  drift, 
we  can  use  the  angular  rate  to  calculate  the  earth's  rotation 
and  carry  out  compensation.  Afterwards,  we  can  find  the 


equivalent  northern  shift  to  be 

£*4-co,,(00)o»at  +  aH,(=«>)sina,*»£^x5sa.-u£ysina,  (12) 

We  then  project  it  on  a  corresponding  axis  and  obtain  the  com¬ 
pensation  quantity  of  the  shift 


Ex  * £*•« ceea, •» £ .cr^a.  —  £  J4n a^osa, 
Ef—  — £'/»sina4*»£tsin,a1-£xcosa.sina. 


(13) 


The  errors  after  we  use  the  above  mentioned  method  to  carry 
out  compensation  are 


&£**  Ex— £*-  -  (£*sina,  4-  £rcosa.)sinci1 

5£r-  £t-£y=*  —  (£ «SiBO.  +  £rC08CtI)c080t, 


(14) 


The  compensated  errors  are  related  to  the  equivalent  western 
drift  which  cannot  be  measured. 

After  the  leveling  of  the  platform,  £?  b=  ^2sinL  is  a 
known  quantity.  When  the  gyroscope's  precision  coefficient  is 
known,  based  on  the  balance  relationship  I(oo)=KG  £2Zb~EZ  °f  the 
control  current  and  azimuth  angle  rate,  we  can  obtain  the  compen¬ 
sation  quantity  of  the  platform's  azimuth  drift. 

BZ-/Cc(2sin  L  -  l  (°°)  (15) 

When  we  determine  the  accurate  position  of  the  runway  or 
other  land  marks  beforehand,  by  using  a  level  viewing  device, 
sighter  or  other  special  optical  systems,  we  can  determine  the 
vehicle's  longitudinal  axis  or  the  azimuth  angle  of  the  platform. 
In  this  way,  we  cannot  only  carry  out  fast  speed  alignment  but 
the  possibility  also  exists  of  calibrating  the  drift  of  the  plat¬ 
form  along  the  horizontal  axis.  In  reality,  when  we  assume  the 
azimuth  angle  measured  and  introduced  by  computer  is  ac ,  from 
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formula  (8)  we  can  then  obtain 


gx»a>X4(oo)-t-  Q  cos  L  cos  a  —  QcosLcosa, 


£Y»oT,(oo)-  QcosLsin  a  -f  QcosLsina, 


(16) 


I 


It  is  very  clear  that  if  there  are  no  errors  in  the  azimuth 
angle,  the  control  value  of  the  balanced  state  is  equal  to  the 
drift  of  the  platform.  When  the  position  error  is  A  a,  we  can 
obtain  the  compensation  error  of  the  drift  quantity  by  using 
the  relationship  of  ac=a+  A  a 


6£x =  Q  cos  L  sinaia 
5£r=“  QcosLcosaia 


(17) 


Naturally,  when  the  precision  of  the  moment  device  is  very  high 
and  the  structural  errors  of  the  platform  are  very  small,  the 
compensation  quantity  is  equal  to  the  drift  of  the  gyroscope. 

It  is  very  obvious  that  the  precision  of  the  platform's  position 
and  drift  compensation  are  mainly  determined  by  the  measuring 
precision  of  the  vehicle's  longitudinal  axis  opposite  the  runway 
or  land  mark  position  as  well  as  the  measurement  precision  of 
the  runway  or  the  position  of  the  land  mark  itself  and  the 
installation  precision  of  the  platform  opposite  the  vehicle.  It 
can  be  seen  from  formula  (17)  that  when  a=0 ,  <S  E=0 ;  when  a=90°, 

4  Ey=0.  This  is  completely  the  same  as  when  carrying  out 
calibration  of  common  platforms.  Further,  the  precision  of  the 
latitudes  introduced  by  the  computer  is  also  related  to  the  pre¬ 
cision  of  compensation.  For  example,  in  overlooking  the  latitude 

errors,  when  L=a=45°,  in  order  to  cause  (5"1  E  and  (S  E  to  not  be 

X  Y 

larger  than  0.02°/hour,  it  is  necessary  to  guarantee  that  azimuth 
error  A  a  ^  12'.  This  data  shows  that  only  if  the  distance  of 
the  runway  or  land  mark  is  not  very  short  and  the  assembling 


precision  of  the  platform  on  the  vehicle  is  also  relatively  high 
will  this  method  have  certain  significance.  Naturally,  when 
we  use  optical  methods  and  are  able  to  cause  the  measurement 
precision  of  the  azimuth  angle  to  reach  about  6',  we  can 
compensate  the  platform's  drift  close  to  0.01°/hour. 


Fig.  5  Block  diagram  of  initial  alignment  loops. 

IV.  Analogue  Computations  of  Navigation  Errors 

In  order  to  obtain  an  even  more  concrete  understanding  of 
the  operating  performance  of  the  azimuth  rate  inertial  naviga¬ 
tion  system,  based  on  the  fundamental  relationship  of  formulas 
(1)  and  (2)  as  well  as  the  corresponding  error  relationship, 
the  calibration  coefficient  error  of  the  gyroscope  and 

accelerometer  is  A  K  =  A  K  =1x10  4;  structural  error  ^.=20"; 

Cj  a  i 

gyroscope  drift  £  -  A  =0 . 01°/hour  and  £  =0.0 2 “/hour;  the 

^  “  _4 

accelerometer's  zero  position  deviation  V  =  V  =1x10  ; 

X  Y  ^ 

initial  velocity  vxo~VYO=(^ '  state  error  Q=  ^YO= 

5x10  4  radian;  the  theoretical  position’s  initial  value 

L  =40°05',  A  =116036',  a  =45°00';  and  when  the  computation 
o  o  o 

position's  initial  value  L  =40°07',  /L  =116°38',  a  =45°03', 

cO  cO  cO 

we  assume  that  the  flight  locus  can  be  described  by  the 
following  parameters: 


The 

6-10. 
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main  computation  results  obtained  are  shown  in  Figs. 


Flight  course. 

Key:  (1)  Degrees;  (2)  Degrees. 


<30 

(2) 


The  time  history  of  position  errors 
Key:  (1)  Degrees;  (2)  Seconds* 


The  time  history  of  attitude  errors 
Key:  (1)  Degrees;  (2)  Seconds. 


The  time  history  of  azimuth  errors 
Key:  (1)  Degrees;  (2)  Seconds. 


Fig.  10  The  time  history  of  felocity  errors. 

Key:  (1)  Meters/second;  (2)  Seconds. 

Fig.  6  shows  the  flight  loci  of  aircraft  taking  off 
from  an  airfield  runway,  spiral  climb,  cruising,  final  return¬ 
ing  and  landing  as  well  as  other  corresponding  calculation 
values.  It  is  very  clear  that  the  errors  between  the  calcula¬ 
tion  and  theoretical  values  change  with  the  different  courses 
and  distances. 

Fig.  7  gives  the  relational  curve  of  the  time  history  of 
position  errors.  We  can  see  that  basically  they  oscillate  using 
the  Schuler  period;  at  the  same  time,  the  range  of  longitudinal 
error  is  larger  than  latitudinal  error  ^L.  This  is  because 

when  increasing  with  the  latitude,  the  longitudinal  value  cor¬ 
responding  to  the  same  type  of  distance  error  must  be  large. 

Fig.  8  gives  the  relational  curve  of  the  time  history  of  the 
platform's  attitude  errors.  It  can  be  clearly  seen  that  because 
of  the  maneuvering  flight  of  takeoff  and  landing,  the  attitude 
errors  of  the  platform  produce  relatively  intense  changes  with 
the  different  azimuth  angles  yet  the  maximum  value  is  still 
within  the  allowable  range.  During  the  Linear  rate  cruising 
stage,  its  law  of  oscillation  coincides  with  the  Schuler  period. 


The  relationship  of  the  time  history  of  the  azimuth  angle 


errors  is  shown  in  Fig.  9.  It  can  be  seen  that  the  changes  of 
the  azimuth  angle  errors  are  quite  slow  and  their  maximum  value 
does  not  exceed  0.06  degrees;  when  we  consider  that  the  initial 
error  is  3',  the  real  change  value  is  even  smaller.  When  the 
latitude  is  relatively  high,  the  change  quantity  of  the  errors 
is  somewhat  larger.  This  also  tallies  with  the  common  laws. 

To  sum  up,  the  output  power  here  must  be  higher  than  when  using 
synchronizers  and  other  electric  machine  transforming  devices. 

This  also  explains  well  the  conclusion  that  the  attitude  errors 
of  the  strap-down  system  must  be  smaller  than  when  using  the 
output  of  the  electric  machine  device  for  the  platform  system. 

Fig.  10  gives  the  relational  curve  of  the  time  history  of 
the  velocity  errors.  We  can  see  from  the  figure  that  during 
takeoff  and  landing,  because  the  velocity  error  is  related  to 
the  spiral  angular  rate,  the  changes  are  relatively  intense. 
However,  calculating  from  the  error  changing  curve  of  the 
cruising  stage,  its  period  also  coincides  with  the  Schuler  period. 

V.  Concluding  Remarks 

Based  on  the  preceding  discussion,  we  can  see  that  the 
azimuth  rate  inertial  navigation  system  is  relatively  suitable 
for  vehicles  with  not  very  large  changes  in  pitch  and  yaw  atti¬ 
tudes.  It  combines  quite  well  the  advantages  of  the  platform 
type  and  strap-down  type  inertial  navigation  systems  and  avoids 
certain  of  their  individual  drawbacks.  This  type  of  platform 
structure  not  only  has  few  components,  is  simple  in  structure, 
light  in  weight,  has  good  reliability  but  moreover  can  provide 
angular  rate  signals  of  the  pilot's  course  damping  channel. 

Thus,  it  omits  a  course  damping  gyroscope  and  it  can  also  give 
course  information  with  relatively  high  precision. 

During  each  alignment,  we  can  carry  out  calibration  and  com¬ 
pensation  of  the  azimuth  drift  of  the  platform.  This  is  another 
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special  advantage  of  this  system.  Moreover,  if  we  use  a  known 
runway  or  land  mark  position  and  draw  support  from  the  level 
viewing  device  or  other  optical  systems,  we  can  also  possibly 
realize  fast  speed  alignment  and  the  calibration  and  compensa¬ 
tion  of  the  horizontal  drift.  This  is  an  area  worthy  of  serious 
attention. 

Because  the  azimuth  angular  rates  of  most  vehicles  are  close 
to  the  angular  rates  of  platform  rotation  in  present  rotating 
azimuth  inertial  navigation  systems,  there  are  no  special  condi¬ 
tions  proposed  for  the  requirements  of  the  system's  computer. 
Naturally,  the  software  system  is  also  close  to  the  common  semi- 
analytical  inertial  navigation  system. 

Finally,  it  should  be  pointed  out  that,  relatively  speaking, 
the  drawbacks  of  the  azimuth  rate  inertial  navigation  system 
are  that  it  is  not  very  suitable  for  fighters  which  perform 
stunt  flying;  because  the  platform  itself  does  not  have  a  degree 
of  freedom  around  the  azimuth  axis,  it  is  thus  not  convenient 
to  carry  out  calibration  and  alignment  of  double  positions  but 
can  only  carry  out  single  position  alignment. 
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OPTIMAL  GUIDANCE  LAWS  FOR  MISSILES  WITH  SECOND  ORDER  LINKS 
Li  Zhongying 

Beijing  Institute  of  Aeronautics  and  Astronautics 
Abstract 


This  paper  studies  the  optimal  guidance  laws  for  missiles 
with  second  order  links  aimed  at  targets.  It  uses  the  prin¬ 
ciple  of  the  minimum  to  separately  research:  1)  when  the  target 
miss  quantity  is  zero,  the  optimal  guidance  laws  of  the 
minimum  control  energy  index;  2)  when  the  final  state  is  a  zero 
control  intercept  curved  surface,  the  optimal  guidance  laws  of 
the  minimum  energy  index.  Finally,  we  obtained  results  similar 
to  those  for  missiles  with  first  order  delay  links. 


Main  Symbols 

X^  relative  position  vector  of  missile  and  target 

X^  relative  velocity  vector  of  missile  and  target 

X-j  relative  acceleration  vector  of  missile  and  target 

X4  rate  of  change  of  relative  acceleration  vector 

u  the  missile  control  vector 

u  the  missile's  optimal  control  vector 

the  missile  '  s  relative  damping  coefficient 
the  missile  body's  intrinsic  frequency 
A’i  conjugate  state  vector  of  state  system 

v,k  Lagrange  multiplier,  undetermined  constant  vector 

tQ,tr  separately  the  opening  and  final  times 

T  time  to  reach  the  zero  control  intercept  curved  surface 

M  time  from  reaching  the  zero  control  intercept  curved 

surface  to  hitting  the  target 

Cd  angular  velocity  of  line  of  sight  between  the  missile 

and  target 

rate  of  change  of  angular  velocity  of  line  of  sight 
(3  real  numbers  selected  according  to  requirements 
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There  has  been  a  great  deal  of  research  [1,3 ,4, 5]  done 
both  domestically  and  abroad  on  the  optimal  guidance  laws  for 
missiles.  A  great  many  beneficial  results  have  been  obtained 
which  have  had  very  good  effects  on  realizing  accurate  guid¬ 
ance.  However,  the  selected  missile  mathematical  models  have 
been  too  simple  and  generally  all  of  them  use  particles  with 
instantaneous  response  for  processing.  Reference  [1]  re¬ 
searched  the  optimal  guidance  laws  of  missiles  with  first  order 
delay  links.  The  results  were  composed  of  proportion  guidance 
with  a  variable  coefficient  and  correction  items  related  to 
rate  of  change  CO  of  the  line  of  sight  angular  velocity  and 
acceleration  a.  However,  the  actual  missile  has  non- 
instantaneous  response  particles  and  is  not  a  first  order  delay 
link-.  The  aim  of  this  paper  is  to  study  the  optimal  guidance 
laws  of  missiles  with  second  order  links.  Because  its  mathe¬ 
matical  model  is  close  to  the  real  thing,  the  obtained  results 
are  beneficial  for  raising  the  guidance  precision  of  the 
missile. 

I.  Formulation  of  the  Problem 

The  relative  motion  dynamics  model  of  the  missile  and  target 
is 


*t**t 

**i# 

**  **J 

*iCO  =*:» 

***** 

**C^o)  * Xil 

<sixs  —  2i<Dxt  +  CDiu 

*«(*«)=  o 

The  performance  index  is 

*Tad1 

rp 


(2) 


Key:  (1)  u"  indicates  the  transition  position. 


--  3 

In  the  formula,  u  fc.  R  and  were  given  in  advance.  We 
assume  the  target  group  is 
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The  present  problem  of  finding  the  optimal  control  laws  of 
dynamic  system  (1)  so  that  the  system  is  guided  from  the  given 
initial  state  to  final  state  (3)  and  moreover  cause  perform¬ 
ance  index  (2)  to  adopt  a  very  small  value. 

Below,  we  will  seek  the  optimal  guidance  laws  for  two  dif¬ 
ferent  types  of  target  groups. 

II.  The  Zero  Miss  Minimum  Energy  Guidance  Law 
The  target  group  is 

X1 (fcf ^ =0 

The  other  parameters  of  the  terminal  are  free. 

Based  on  the  principle  of  the  maximum  [2],  the  problem's 
Hamiltonian  function  H  is : 

H  =  -g«r  u-i-xr' +  — +  co*u)  (  4  ) 

M 

In  the  formula,  A^,  '  ^3  an^  A4  are  the  conjugate  state 

vectors  the  dynamic  system  (1). 

f  A,(  O-  0 

AS(0-“V(0  (5) 

Aj(  t  )—-/;(  •*  )-W:X/  /  ) 

■  A.(  t  )  —  -Aji  :  )  -2iU>X.C  <  ) 


Their  cross  section  conditions  are 


(6) 


I  0 

j  a.j( *») =  o 

L  X4(*,)-  0 


Based  on  the  principle  of  the  minimum,  the  optimal  control  laws 
and  corresponding  conjugate  state  parameters  should  cause  the 
H  function  to  take  the  minimum  value  and  from 


aH 

Mi 


=  0 


we  obtain 


i(  =  (7)' 

We  can  see  from  formula  (7)  that  to  determine  the  optimal 
control  laws,  we  must  determine  conjugate  state  vector  -4^  ( t) . 

From  the  first  two  equations  in  equation  group  (5)  and  the 
corresponding  boundary  conditions,  we  can  solve 

X.CO-* 

x,(  o- v  a-  * ) 

Therefore,  we  have 

’  Aj(  <)  =  -  v  (<,-  t  )+®*X4(  O  (8) 

We  then  seek  a  derivative  of  the  fourth  equation  of  equation 
group  (5)  and  obtain 

a4(  i  )=*-*,(  i  )+26®a4(  «  ) 

,sr.  we  substitute  in  formula  (8)  we  obtain 


A.(  0-24®A«(  t  )  +  <*%(  t  )-  V  (/,-  t  ) 


(9) 


We  can  know  from  formula  (7)  that  to  find  the  optimal  control 
laws  it  is  only  necessary  to  solve  ^4  (t)  .  In  equation  (9), 
we  let 


and  therefore  equation  (9)  becomes 


A«(  *  )  ?  )  +  <a1X4(  t)  =  vt 


(9)' 


When  =0,  t=tf  has  /t'|0)=0  and  moreover  we  suppose 
^t.(0)=0.  Formula  (9) 1  is  a  second  order  non-homogeneous 
linear  differential  equation  with  a  constant  coefficient.  We 
can  find  its  solution  by  using  the  Laplace  transformation: 

i 

i 


K  (  * )  - 


VT 

CO2 


_2|V 

Cl)3 


— ^r-[2|co8flT  +  —  ^.^.t-sinQTj 


(10) 


In  the  formula, 


Q-fflO  -VV" 

By  substituting  formula  (10)  into  formula  (7) ,  we  can  obtain 
the  optimal  control  law: 

i(  D  —  v  {(*,-  t  )-- g-  +  [2$cosQ(f,-  <  ) 

+(12-7t)‘7isiTiQ(f/~  °]}  (ll) 

For  the  final  determination  of  the  optimal  control  law,  we 
must  determine  constant  vector  v  and  for  this  reason  we  must 
substitute  formula  (11)  into  the  system's  state  equation  group 
(1)  to  solve  the  state  parameters. 

From  the  third  equation  in  equation  group  (1)  we  have 

x3 ( t) =x4 (t) 
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When  we  substitute  the  fourth  equation  in  equation  group  (1) 
into  the  above  equation,  we  obtain 


*j(  l  )- -©**,(  t  )-2l®*t(  /  )  +  co15(  t  ) 

We  then  substitute  the  third  equation  in  equation  group  (1) 
into  the  above  equation  and  obtain 

*s(  t  )  +  2|a»,(  t  )+«**,(  t  )-c o*u(  t  ) 

We  substitute  formula  (11)  into  this  equation  and  obtain 

*j-r240J*,  +  a),Xj  *  i  ) — +  ^24cosQ(f,-  t  ) 

+  (  Q  Ot-  t  )  ]}  (12) 

This  is  also  a  non-homogeneous  linear  differential  equation 
with  a  second  order  constant  coefficient  which  can  still  be 
solved  by  the  Laplace  transformation.  Finally,  we  obtain 

*,(/)--  v  <(ff—  «)-[*>  +  ]  e‘"'siii (£#  +  *,) 

~4<D(At'jT7?e't*<,r>3inCQ(f^ 

+^y~ipe*‘-<'/*‘>c°sCQ t  )  +  +,)}+*„  T^re*t-'cos(Qf  +  ^)  (13) 


In  the  formula 


■ii  ”  tg‘l(  |  ^|rtg(Q// + w] 

*4{l.‘(i^lP"(°l,**J 


a  — 


b  ~  44a>*(  1  -V)  ,/*sin  +  ^ 


They  are  ail  the  missile  body's  known  characteristic  parameters 
%y  and  CC>  as  well  as  the  function  terminal  time  tg.  When  t^ 
is  given,  it  is  an  unchanging  constant. 


Based  on  the  first  equation  in  equation  group  (1)  and  its 
initial  conditions,  by  substituting  in  formula  (13)  and  inte¬ 
grating  we  can  obtain 


*t(  t  ^ C«'l",acoa(Of«  +  \|i, 4- -e'^'coaCQf  +  *|i.+  40) 

-  ’  {(«  -  f )  H *  (  « fl ' *  « 

+  k  *  }  g-w'^cos(Qf  +  ^  +  ^)+  4<afc  ■  _%ty  rsinCQ (*,-  t  )+4>,~4>J 

-I Oami<Bttrn 
4CD*(  1  — Q  (*/*”  *  )  ^0  +  ^*3  ^  — 4i)1  lC0S^^  (Jt~~  O 

+  *.+W-(<A— ^-)+C,«h#,f  4,cd)1  (14) 


In  the  formula 


CXt"U,  +  j— ]'  *  (Qf,  + 1|>,  +  •>!>«) 

_  ^±pH^,r.,>cosiQia  +  ^  +  _-_L_<,-t.<.,«,)sinCQ(/,-/,) 

O  p'ivl,r‘9> 

+  *.  -  +  J  -  C  Q  0, -  O  4-  M>«  +  +4  J  >  “  4t^3(  i  irpjTTtCosC  Q  0,  -  /,) 

+  $.  +  **] 

w(£] 


Based  on  the  first  equation  in  equation  group  (1)  and  its 
initial  conditions,  we  substitute  in  formula  (14),  integrate 
once  ^nd  obtain 


*|(  O -*!.  +  *„(  <  +  )[«•*“«(  t  -f,)a»(Q#,  +  *t*t>4) 

- ^ — COs(Qf04-Hi04-2^j4--^-cos(QY+iJis  +  2\i)J  j 

r  ’  k  *4--4o-i**i^n  *?*<*+*■+«« 

(a1  +  bz)wt  -»•<»■*•> 

-  ^  CH  -4»  ih.  4>  ih  h  - f 1 _ _  , 


/  ^  W  ^ 

”co5(Q/  +  ^,  +  i|i«)+-j2 S(Y~4iyrTSintQ(/f-  <  )  +  %-  2iO 
~4^(l-l;)‘ -siT,CQ(f>~  ‘  ^  +  2^3 + ^g4( t , t)l  ;cosC  Q  0,-  /  )  +  ^,  +  2-W 


r 4 


In  the  formula 


c,a,<„  i.oo=.-c,(m„  4,co)f0+  -*£ 

(a1  +  b1V* 

x  sin  CQfe  +  tb,  +  2vJ  +  ^  -r4-V - **»-<*rVcoi(Q*t  +  *,  -Hji.) 


oo- 


-  4m»(1~  4»)  ■  :sin  CQC#f  —  O  +  -  2^4)  +  J^(JITp)T7-isin  CQ  (</  —  /*)  -4-  il>e  4-  2U-J 

Q  .l.(,r,0, 

-  4Ta«("i  -  i») 1  »cosC  Q  C»> - O  +  *,  +  2^,3 
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From  we  can  determine  constant  vector 


v , 


;  [  *l»  +  (f / -  *,)  +  [<•*•*»(< » -O003(Qf,  +  tj),  +  *j»4) 

- — cos(Q*0  +  ia  +  2i»4)  +  cos(Q/,  +  *.  +  2t*>*) j  x  {  —  3  °— 

”  [  <*  +  ~"a>^-~]  ^^rsin(Q<i+ +i  +  2i(iJ  -  - — ■e'u“',cos(Q/.  +  +  t|>4) 

1  O 

+  4<a3(  1  -il)‘  *®n ~4<B*(l-lt)l'lSin(li’,  +  2,i,‘)4'4|'o,C  1  —  4* ) 1  s 


xcos(<j»,  +  2^4)+C,(/,, /„  4,  c o)/,  +  C, (/„/,, 


(16) 


V7e  substitute  the  v  value  into  formula  (11)  and  obtain 


i(  t  )*  —  |x10  +  jt.0  (<,-/,)  +  ^  [g'***^^  - <.)cOS(Q/0  -M|l0-r  TtJ 

- - — cos  (Q;„  +  Al>t  +  2^4)  +  — - — cos(QY,+t*,-i-2Tj)4)|  x  \(1,~  t  ) 

, )  +  « )j|x {2^21 

-('!+  nt»‘r  ^5i”  «»i+ ♦> + *0 +*,**.) 

g-U“‘/  J  Q 

•  +4o»(1_^)t'«  SiD (*,  +  2lil.)  +  4ioi«(  J  _;l);  » 


xa»(4,  +  2^J  +  Ci(ihU,  i,  cD)/f+C,(fft/M  £,u>) 


)•* 

i 


(17) 


Letting 


1 


X  e'u*‘t<xa(Qtf  +  +  *4)  +  x  4 tsin ( *,  -  2ij>4) -  4<a»(  L  ZTipr 

x  sin (ifr,  +  2*4)  +  4fo7( Y " |>) i/iCwC^  +  +  C, (** l,  w)*t 

+  C.(W«  CD) 


When  tf  is  given,  F(tf/tQ,  %>  ,  t£>  )  is  a  constant  quantity 
determined  by  the  missile  body's  parameters  4>  /  and  time  tQ. 
In  this  way,  formula  (17)  can  be  simplified  into 
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I(  t  )=»  - 1  a  ^  ('■)  [*“'*(<> +  tl>4) 

— «08(Q#.  +  t>t  +  2^t)+-^j-^-co8(Q#f+  ♦i  +  2't«)|  x  (/,  —  t  ) 


It  -l •<',-»)  c  Oi1—  1  1 ) 

-ir+ [2SotQ  o+ir-n"*™-  ‘ >]} 

»j(l,-|')'tf(l,1.  5. ■»)}"' 


When  t=>tQ,  optimal  guidance  law  is 


i(0  »  ~  {  <„  +  *».(*,-*.)  +  ^ (  ft?)  [«'*”*(*>  -  Ocm(Q  U  +  *•+  40 

- — — cns(Qf0  +  ij>(l  +  2ii)4)4--^ — cos(Qf,  +  ifc,  +  2TiO|  *  'fit-* o) 


24 

CO 

®  l 

24*-  1 


1  v  K<>-0* 
)'  1  3 


+  F(f„/.,  4,®)} 


We  let 


K I ,  « )  -  ^_  +  J  jVi)*n  CQ(*,-0+W  } 

x{-L~f—  +  FO"  U,  4,®)}“ 

i  r  »-*•*» 

I  </„/„  4,  + 

x cos (Qf,  +  <fr,  +  2^4) H — ^—cos(Q»,  +  ’•{),  —  240  ]  x  - ~ 

+  m(  r 3ir)3ii»CQ(/r-0  +  t|,,3  { x  { +F(i„U,  4,®)  • 

f 

In  this  way,  formula  (19)  is  simplified  into 

4,  ®)C*.»+x., (#,-*,»-  £(**#..  4.®)*,,  (20) 

The  guidance  law  shown  by  formula  (20)  can  use  the  appro¬ 
priate  selection  time  t^-t  which  can  change  it  into  proportional 
guidance  with  a  variable  coefficient  (in  proportion  to  the  line 
of  vision  angular  velocity)  and  the  guidance  law  [1]  composed  of 


correction  items  with  one  item  related  to  the  acceleration  and 
line  of  vision  angular  acceleration. 


When  CD  ->  oo,  the  dynamic  response  process  of  the  missile 
is  lost,  that  is,  the  response  of  the  system  is  without  oscilla¬ 
tion  as  well  as  without  delay.  At  this  time,  F(t,,t  ,  £>  , &  ) -Q . 

to 


K  ( it,  itf  4  .  m  )  - -jr  _3f  jT 

Ui„U,  4,  ®)-  0 


Therefore 

v_  3  C^ie  (.^t  C?]') 

This  is  the  optimal  guidance  law  of  the  instantaneous  response 
particles.  After  selecting  the  different  t^-t  values,  we  can 
obtain  the  proportional  guidance  law  of  the  various  different 
systems  [3 ] . 

III.  Guidance  Laws  of  the  Target  Group  as  a  Zero  Control 
Intercept  Curved  Surface 


When  dynamic  system  (1) 
intercept  curved  surface  L, 


is  guided  towards  zero  control 
its  target  groups  is  [1,2]: 


j  jr.cn+nx.cr)- o,  n>o 
1  jc,(r)=*  o 


(22) 


It  is  assumed  that  the  introduced  time  T  of  the  zero  control 
intercept  curved  surface  L  is  already  selected;  time  JU.  from 
entering  the  zero  control  intercept  curved  surface  to  hitting  the 
target  is  undetermined. 


The  main  reason  why  x^{T)=Q  is  necessary  is  because  it  is 
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required  that  we  cause  the  system  to  be  able  to  be  maintained 
within  this  curved  surface  after  bringing  the  system  into  the 
zero  control  intercept  curved  surface. 


The  performance  index  is  still  formula  (2)  and  at  this  time 
Hamiltonian  function  H  of  the  system  does  not  change.  Thus, 
the  conjugate  state  set  of  equations  of  dynamic  system  (1)  is 
still  described  by  formula  (5) .  However,  the  cross  section 
conditions  change  to 


X,(7)  =  v 

>.t(r)=v4 

x,(D-  k 
■  x«(D-  o 


In  the  formula,  v  and  k  are  Lagrange  multipliers  which  are 
undetermined  constant  vectors.  Based  on  these  boundary  condi¬ 
tions  solving  conjugate  state  set  of  equations  (5) ,  we  obtain: 

X.(  t  )-  v 

x»(o-»(r  +  i»  -  t ) 

A,(  /  )--  *(7  +  t  )+<B1X4(  i  ) 

A«(  t  )=*2i<BA4(  t  (  0+  *  (7  +  V-  -  t  ). 

Based  on  the  results  in  Section  n,  we  obtain 


Ul)-yjtr+»-l)- -f- 

2il—  1  T) 

+  "(  I’rjijTTJ-sfoSH  7  +  ^  —  I  )J  j- 


[24cosQ(7+  I*  -  t  ) 


When  we  substitute  /^(t)  into  formula  (7),  we  obtain  the 
optimal  guidance  law 

“(<)--  v  •[ (7A  n  -  O — ^-  +  ^’i*,T**w>[2£cosQ(7+  4  -  /  ) 

+  t  .’"'iiVTT'  sin  Q  (  7  +  V-  -  t  )  V:  (2S) 


When  we  substitute  formula  (25)  into  dynamic  system  (1) 
just  as  in  Section  II  we  can  solve  x^(t),  x^(t)  and  x^(t). 

{(  T  +  *  -  t  )  -  [(  T  +  *  y +  ^  7  y  )4~ 1 )  1 1  ] ' ‘  V*-»n(Qf  -i-  ♦,) 

~4^(i-&iF»ginCQ(r+tt  "  *  )-W-(ai  +  6I)‘  V‘--(^'> 
xcos(Q*-Hk)-  4ffl(.  -  j.|>)1/>3inCQCr  +  ^  -  *  )  +  <!>, 3 


+  4rff(- 1  r^-.-cosC Q  ( r  +  n  -  n+wJ+y^^'cMCOf-r*,)  (26) 

jct(  t  )=xit+  [g't"**a»(Q7.  +  ii>t-i-t|i4)-g'l"<cos(Q<+^,-)-t|)i) 

-  v{[(r^-^)  f 

(er  +  h1}1'1  -  i.<r»»*o 

+  (— ^  ~  g--—  cc<Q* +4.+ *.)  +  isrcr-^ 

xsinCQCT+f1-  t  )-rtfr,-*J — 4m.(.  -  _%iy  ..tSinCQCr  +  ^  -  <)+$,  +  $  J 

f  fl  -I 

+  4ff(1  "■  -^ticmC Q(T+  ^  -  o  +  +,  +  *j-  j/r  +  li)*.-  ~r. 

+Ct(r  t  f1,/,,  4.  00 )j  (27) 

x,(  t  )3*x„4-jfl0(  t  ~ij+  f  -#.)«*(Q*. +  +•  +  *«) 

- - — cos(Q<#  +  ifr,  +  2i>  J  +  -^-a»(Qf  -r  ’4>0 1 2’K)  J 

,  v  |[(r + ,)|.|].[(r+ »).+  a^xy  )±i)f]-'  ^sin(Q<^^2W 

xsinCQCT  +  H  -  /  )  +  *,-2^«D  - 4~ , ^j-^-|iynsiii CQ(T  +  ^  -  <  ) -i- ij>0  +  2ij>  J 

+4Sff( ^r^pncosCQCr  +  i*-  0  +  4>,t2’W  -[  (r  +  n)f, — f-J  < 

+0,(7*  + 4,®)<  +U±!1.>fl  -  -if-+C,(  T+P-J,,  6,®)}  (28) 


J  ,  x. 


From  formula  (28) ,  letting  t=T  we  can  obtain 

*i(r)-x?(D- v  {( LrW  +  +^T,  n,  *,->){•  (29) 


In  the  formula 


*!(  T) 


- L- - COS(Qf0>— - —  CO$(QT  f  <j)<,-r2liOj 


IX AT,  S,  .  0> )» - [( T  +  H  )•+  11  ^*n(Qr+*,  +  2*J 


~  (~fll ) 1/1  e'‘*(1T*“)cos(QT 4-  +,  + +4)+ 


i«C2T*“> 


4co3  (  1  -i2)1-1 


gl— 

Xsin(Q^  +  iii„-2Tij4)  — 1  ^v'tyrt-siTi  (QP-  4-  4,  -r  2tji4) 

+  4-4— ^  —  T^ncosCQM- + 4>,  -  24,)  ~C.  C  T  T  I*  ,  l,co)T 

-*-c0(r  +  , f8, 4 ,  co ) 


From  formula  (27) ,  letting  t=T,  we  can  obtain 

*t(r)-x?(n-  v  H(f  -t,)  +  *Xt(T,  H,  I,  C0)1  (30) 

In  the  formula 


*8(7*  )  -*M  +  --(^a_|t)Cg-i*,’cos(Qf,  +  4.  +  4J  -  e*l-Tcos(Q7*+  4,  -  v.) ) 

AJf.(7\n,|,<p;«[(r+  »)l+  +  1->' j‘  '  ^siT>(Qr+  4.  +  4J 

+  -(g-^  ~  g't"<,T*“,cos(Q7’+  4,  +  *4)+ 

e-i« 

xshiCQix4-4ii~4J  —  4CQ2(-"1  sin(Qt*-f-  4»-i-4«) 


+  aci/(‘ \  ~ it)T;jcos(GH  +  ->• 4J  +  C, ( T  +  H I,  «) 


Because  the  system  should  reach  the  target  group  when  t=T 
it  should  satisfy 

xl(T)+  /Ux2(T)=0 

From  Formulas  (29)  and  (30) ,  we  can  obtain 


*t(  n +•»*}(  r )- v 


f  (T-UY 
(  3 


■+  ncr-o+^cr-oti.^cr,  n,  i, ©) 


+  V-&X.XT,  4,  <*>)}  =  o 


We  can  find  v  from  this  formula. 


( T  -;t)3 


_ *f(r)+4*°(r) _ 

*(T -1Jl+vl(T-u)  +  bXXT,»,i.  *)  +  *lXXT,  n,  4.®) 


When  we  substitute  formula  (31)  into  formula  v_i; . 
obtain  the  optimal  control  law: 


u(0--c*I(7’KM(r))|(r  +  n-  o-  + 


^cosqct-1-  m.  -  t )+(  ^^.^i-zSinQc r -  m.  -  t  )jj  •  | 


",  (  (T-uY 

J )  (  3 


+  ^  (7*  +  A.V,(7\  n,  4,  a>)~*lX.(T,  v-,  4. 


When  t=tQ,  the  optical  guidance  law  is 

•CO  -  -  c*j(  r) +•»*}( r)]{  c  r  +  n  + ^7  1 

xsiT)CQ(r  +  ^  -<,)  +*,3}  +n8(r-o 


+d^t(r,  n,  i.<o)+^^,cr,  ti,  1,  cd) 


The  above  formula  can  be  rewritten  as 


-kxt ,  4,  cD)cx„+x:o(r  +  m. 

-*,)!! -£<.(r,  H,  4,  co)x„ 


In  the  formula 


-  *■  «'■  /fch '** »-**». '**.*  .*■  .  *.•  ,•  ’■*.  V  '  1.  r _  -  *..  ‘0  V\*Va‘.-.', 
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Kl(T,  l ,  O )  =  {  (  T  +  i*  -  O - -f-  +  -(~t Ly} 

x  sin  (  Q  (  7  -r  (*  -#,)  +  +,]}  {  +  >*  (  7  -  0 1  +  ^ (  7  -  /„) 

+  A*,(7,  (*,  i,  a>)-r!iA^.(r,  I1,  I,  ®)}  ' 

Ll(.T,  I*.  &.  «>)-— («' '*‘*(7  +  E  -Oc«(Q*.  *+.  +  +>) 

—  ^e'i",cos(Qi,  -t-  4>o  +  2+ .)  +  -^■<?'>Tcos(Q7 -r  +.  +  2+J 


-(**•‘•'0)8(07-.-+.++.)}  (  (-T-zuy  T(*(r-f.),+t*icr-f,) 

+aa\(7,  i*.  i,  co)+(iAX(r,  (*,  4 ,  co ) } 

The  "^2.'  "^3»  a'  5:3  an<^  t^ie  tf  section  should  all  be 

replaced  by  T  tpo. 

In  optimal  control  law  formula  (33)  or  (34),  the  is 
still  undetermined.  Only  when  the  is  determined  can  we  cal¬ 
culate  and  completely  determine  the  optimal  control  law.  For 
this  reason,  we  used  the  condition  of  x^(T)=0.  Although  this 
condition  has  already  been  considered  in  the  limit  value  con¬ 
ditions,  yet  because  it  is  established  in  the  optimal  control 
process,  this  optimal  control  law  is  unrelated  to  A-  and  there¬ 
fore  the  condition  of  x^(T)=0  is  still  not  really  used.  In 
order  to  guarantee  not  separating  from  this  curved  surface  again 
after  introducing  the  system  into  the  zero  control  intercept 
curved  surface,  we  must  select  the  JUs  in  this  way  so  that 
x^(T)=0.  From  formula  (26),  letting  t=T,  we  can  obtain: 


By  substituting  in  formula  (31),  we  can  obtain 

-cxtcn+njcscnj  -[cr ^ *y+  (^°(r ^ )~1)~j  e'^sincor+v,) 

.  1  (II  (  *  J"  +  U  ) 

-  -  /;  \-;Tl  ,S1T»  (Qn  +  *.)  -  (a1  +  6:) 1  ■  V*-‘«^,cob(Q7-  +  V,) 

4co(^  i  — -  ; 

-*•*  Oe~iatl  N ) 

- iSfi^FF15" (S“ **°+ +*•>> 

x  +  u<r-i,)*^i‘,(r-(.)^iJf1(r,  >■,  s, ») 

1.  0 

+  HAX.(r,  n,  I,  +T^r^“Tcos(Qr+^)«  o  ( 


This  is  a  high  order  transcendental  equation  which  can  use  a 
numerical  sol  fion  to  solve  yU-.  There  are  possibly  many 
values  to  satisfy  equation  (35).  At  this  time,  we  should  take 
the  smallest  one  among  the  M'  which  is  larger  than  zero. 
Naturally,  it  is  also  possible  that  no  exists  which  can 
satisfy  equation  (35).  Under  these  conditions,  the  system  is 
guided  to  the  zero  control  ir.tercept  curved  surface,  yet  it 
cannot  cause  the  system  to  be  maintained  in  this  curved  surface. 
At  this  time,  it  loses  significance  for  the  guidance  of  the  zero 
control  intercept  curved  surface.  However,  physically,  control 
to  cause  the  relative  acceleration  be  zero  exists  and  therefore 
the  equation  should  have  a  solution. 

When  £)  —too,  the  characteristic  loss  of  the  missile's 
second  order  link  becomes  particles  with  instantaneous  response. 
At  this  time,  (T,^,  %,£>)  and  L  (T,^,^,,  to)  of  formula  (34) 

become 
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T  M.  -u _ 

KL(T  ,  ^  *  4,^0"“  /  'T'  _  J  'J  5 


U(T,  n,  I,  “)=  0 


Optimal  guidance  law  u(t  )  changes  to 


i»(T  +  4  -.',) 


*-o 


(36) 


This  is  the  result  obtained  by  guiding  the  instantaneous  re¬ 
sponse  particles  toward  zero  control  intercept  curved  surface 
L  [3]. 


IV.  Concluding  Remarks 

Above,  we  obtained  the  optimal  guidance  laws  of  formulas 
(20)  and  (34)  for  missiles  with  second  order  characteristics 
for  two  types  of  target  groups.  These  two  formulas  are  sim¬ 
ilar  in  form.  If  the  missile  1  s  acceleration  is  expressed  by 
the  rate  of  change  of  the  line  of  sight  rotation  angular  vel¬ 
ocity  and  angular  velocity,  and  we  let 


T  -r  n  -U 


^•i  o  (G*m  ~  3*;c) 

4c  0JC,Q) 


then  formulas  (20)  and  (34)  can  be  expressed  as  [1]: 


X.j,  ®  *  0*  4  »  ®  ^  ) 

x..,,,  a,  0  ,  4,  CO,  M>)- 


cf,Cajc.<-rir.a) 


x  (gjC|^ *„ -i-fi, x (ax„ -f- 3^,)j -r [AT^x.a, a,  0,  |  o,  a) 

X  xt^x~^M~2Ll(-X'°'  *m’  “*  3’  **  ®'  ^ 
x  *{,(«!*+ 3*„)  fo.;  ( 


In  the  formula,  and  CO  are  separately  the  angular  velocity 
of  the  line  of  sight  and  the  angular  acceleration;  a  and  /$ 
are  real  numbers  and  can  be  selected  according  to  requirements. 

It  can  be  seen  from  formula  (37)  that  the  optimal  guidance 
law  for  missiles  with  second  order  characteristics  is  composed 
of  proportional  guidance  with  a  variable  coefficient  and  a 
correction  item  related  to  the  acceleration  and  line  of  sight 
angular  acceleration  under  specified  transit  time  conditions. 
Naturally,  to  realize  this  type  of  guidance  law,  aside  from 
needing  to  measure  the  angular  velocity  of  the  line  of  sight,  it 
is  also  necessary  to  measure  the  rate  of  change  and  acceleration 
of  the  line  of  sight  angular  velocity.  Although  the  computation 
of  its  coefficient  is  relatively  complex,  yet  they  are  already 
known  functions  of  -x^g,  x2o'  a'  $  >  ,  Co  and  .  We  believe 

they  are  not  difficult  to  complete  in  today's  flourishing  of 
electronic  technology.  Because  the  mathematical  models  selected 
for  the  missiles  went  a  step  further  than  the  mathematical 
models  provided  by  the  listed  references,  we  will  thus  have  even 
higher  guidance  precision. 
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